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The analytic quadrupole octupole axially symmetric model, which had successfully predicted 
^^®Ra and ^2®Th as lying at the border between the regions of octupole deformation and octupole 
vibrations in the light actinides using an infinite well potential (AQOA-IW), is made applicable to a 
wider region of nuclei exhibiting octupole deformation, through the use of a Davidson potential, + 

(AQOA-D). Analytic expressions for energy spectra and B(E1), B(E2), B(E3) transition rates 
are derived. The spectra of and 224,226rpj^ described in terms of the two parameters 

00 (expressing the relative amount of octupole vs. quadrupole deformation) and do (the position 
of the minimum of the Davidson potential), while the recently determined B(EL) transition rates 
of ^^'^Ra, presenting stable octupole deformation, are successfully reproduced. A procedure for 
gradually determining the parameters appearing in the B(EL) transitions from a minimum set of 
data, thus increasing the predictive power of the model, is outlined. 

PACS numbers: 21.60.Ev, 21.60.Fw, 21.10.Re, 23.20.Js 


I. INTRODUCTION 

Rotational nuclear spectra have long been attributed 
to quadrupole deformations [T]. However, octupole 
deformations [corresponding to reflection asymmetric 
(pearlike) shapes] [5H1] are supposed to occur in cer¬ 
tain regions, most notably in the light actinides [5HH] 
and in some light rare earths HE]. The hallmark of oc¬ 
tupole deformation is a negative parity band with levels 
L'^ = 1~, 3“, 5“, ..., lying close to the ground state 
band and forming with it a single band with = O’*", 
1“, 2+, 3“, 4+, ..., while a negative parity band lying 
systematically higher than the ground state band is a 
footprint of octupole vibrations. 

The transition from the regime of octupole vibra¬ 
tions into the region of octupole deformation has been 
considered by several authors [HHHj- In the analytic 
quadrupole octupole axially symmetric (AQOA) model 
the actinides lying on the border between the re¬ 
gions of octupole deformation and octupole vibrations 
have been described, making the following assumptions. 

1) Quadrupole and octupole deformations are taken 
into account on equal footing, their relative presence de¬ 
scribed by the only free parameter in the model, 0o- 

2) Axial symmetry is assumed, in order to keep the 
problem tractable. 

3) Separation of variables is achieved in a way analo¬ 
gous to the one used in the framework of the X(5) model 
m. describing the first order shape phase transition be¬ 
tween spherical and quadrupole deformed shapes m- 

4) An infinite well potential is assumed appropriate 
for the description of the border region, as in the E(5) 
[m and X(5) [T6| models, the former one describing the 
second order shape phase transition between spherical 
and 7 -unstable nuclei. Therefore we are going to call 


this solution the AQOA-IW model. 

A different approach to the problem of phase transition 
in the octupole mode has been developed by Bizzeti and 
Bizzeti-Sona muni, characterized by the introduction 
of a new parametrization of the quadrupole and octupole 
degrees of freedom, using as intrinsic frame of reference 
the principal axes of the overall tensor of inertia, as re¬ 
sulting from the combined quadrupole and octupole de¬ 
formation. The main differences between the two models 
are: 

1) The AQOA model is analytic, while the model of 
Refs. [E |20] is not. 

2) In the AQOA model the symmetry axes of the 
quadrupole and octupole deformations are taken to coin¬ 
cide, in order to guarantee axial symmetry, while in the 
more general framework of Refs. [EllSI nonaxial con¬ 
tributions, small but not frozen to zero, are taken into 
account. 

In both models da mi HO], ^^®Ra and ^^®Th appear to 
lie close to the point of transition between octupole de¬ 
formation and octupole vibrations, with heavier isotopes 
corresponding to octupole vibrations and lighter isotopes 
exhibiting octupole deformation. 

The recent experimental verification of stable octupole 
deformation in ^^'^Ra [5T| stirred interest in octupole de¬ 
formation in the light actinides and their theoretical in¬ 
terpretation. The AQOA model can be made applicable 
to deformed nuclei near the transition point by replacing 
the infinite well potential by the Davidson potential [22] 
of the form /3^ -|-/3 q// 3^, which contains an additional free 
parameter, the position /3q of the minimum of the po¬ 
tential well. The flexibility acquired through the replace¬ 
ment of the infinite well potential by the Davidson poten¬ 
tial has been demonstrated and exploited in the case of 
quadrupole deformation in |23j . The analytic quadrupole 
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octupole axially symmetric model with a Davidson po¬ 
tential, to be called the AQOA-D model, is the subject of 
the present work. In addition to the spectra of 
[3 and 224.226rp]^ puns], the recently measured PT] 
electric transition probabilities of ^^^Ra provide an excel¬ 
lent test ground for the model, already exploited in the 
Bizzeti and Bizzeti-Sona approach |26j . 

The above mentioned work on the octupole degree of 
freedom has been developed in the framework of the col¬ 
lective model [T] . Alternative approaches include the fol¬ 
lowing. 

1 ) A complete algebraic classification of the states oc¬ 
curring in the simultaneous presence of the quadrupole 
and octupole degrees of freedom has been provided in 
terms of the spdf-interacting boson model PTHSO] , which 
has been successfully applied to Ra El, Th PH, u pa, 
and Pu pa isotopes. Mean field studies of the critical 
point for the onset of octupole deformation in quadrupole 
deformed systems have been carried out in Refs. P51IM] . 


II. THE ANALYTIC QUADRUPOLE 

OCTUPOLE AXIALLY SYMMETRIC (AQOA) 
MODEL 

A. Formulation 

In the AQOA model [15] the following assumptions are 
made: 

a) The axes of the quadrupole and octupole deforma¬ 
tions are taken to coincide. In other words, axial symme¬ 
try is assumed, while the 7 degree of freedom is ignored. 

b) Levels with K ^ 0 (where K is the projection of 
the angular momentum on the body-fixed z' axis) are 
ignored, since they are lying infinitely high in energy PU- 

The Hamiltonian of the AQOA model reads 

H = - V ^ ^ ^3 ^ I 

2Bx PI dPx^^dPx (S{B2PI + 2R3/3|) 

+ m,/33) (1) 


2 ) An alternative interpretation of the low-lying nega¬ 
tive parity states in the light actinides has been provided 
in terms of clustering [85H39j . The recent experimental 
findings for ^^"‘Ra [H] seem to point against this inter¬ 
pretation, but wider evidence in more nuclei is desirable. 


where P 2 and P 3 are the quadrupole and octupole defor¬ 
mations, B 2 , B^ are the mass parameters, and L is the 
angular momentum operator in the intrinsic frame, taken 
along the principal axes of inertia. 

The solutions of the Schrodinger equation read pa 


3 ) Relativistic mean field calculations involving the oc¬ 
tupole degree of freedom have been carried out both in 
the light actinides region [401 ST] and in the light rare 
earths PUSS], corroborating imsi] the transition from 
octupole deformation to octupole vibrations in the light 
Th isotopes. 

4) A hybrid approach combining the algebraic ap¬ 
proach of the interacting boson model of 1 ) with the rel¬ 
ativistic energy density functional theory of 3) has been 
recently developed jH] and applied the Ra and Th iso¬ 
topes pu SH] , as well as to the rare earths Ba and Sm 
[4bj , again corroborating the transition from octupole de¬ 
formation to octupole vibrations in the light Ra and Th 
isotopes. 

5) The extended coherent state model (ECSM) has 
been successfully applied to the description of several 
negative parity bands in Rn P7], Ra PTHSU] . Th [50] . 
U [50], and Pu [50] isotopes. 

The AQOA-D model is described in Section II, while 
in Section III numerical results are provided and subse¬ 
quently discussed in Section IV. The integrals needed in 
the calculation of electric transition probabilities are cal¬ 
culated in Appendices A-C, while in Appendix D some 
details of the derivation of the Hamiltonian, the method 
of solution and the comparison to other approaches are 
given. 


<J>!(/?2,/33,0) = iP2P3)-^^^^tip2,P3)\LM0,±), (2) 

where 9 are the Euler angles describing the orientation of 
the body-fixed axes x', y', z' relative to the laboratory- 
fixed axes X, y, z, while the function |LM 0 , ±) describes 
the rotation of an axially symmetric nucleus with angular 
momentum projection M onto the laboratory-fixed z-axis 
and projection K = 0 onto the body-fixed z'-axis [T] 


|LM0,±) = ± (-I)^)DoV(^)> 


(3) 


with I?(0) denoting Wigner functions of the Euler angles. 

Wave functions with the -I- label correspond to positive 
parity states with L — 0, 2, 4, , while these with the 

— label correspond to negative parity states with L = 1 , 

3, 5, .... 

The Schrodinger equation can be simplified by intro¬ 
ducing [5n 15^ 


P2 




B = 


B2 + B3 
2 


(4) 


reduced energies e = {2B/h?)E and reduced potentials 
V = {2B/fP)V, as well as polar coordinates (with 0 < 
P < 00 and — 7 r /2 < p < 7 r/ 2 ) [511 [5^ 


P2=pcosp, P 3 = Psmp, p=\lpl + pl, ( 5 ) 
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leading to 




1 ^ 


+v{j3,(j)) 


L(L + 1) 
3;9^(1 + sin^ 4 >) 
3 

sin^ 2(/) 


Cl 


1 52 
^32 ,9(^2 

vE-±(;9,<^)=0. (6) 


In addition, separation of variables can be achieved 
by assuming the potential to be of the form w(/3, (j)) = 
u{/3) + where is supposed to be of the 

form of two very steep harmonic oscillators centered at 
the values ±(j)o- In this way Eq. ([^ is separated into 


r 92 

1 d 1 

/ L(L + 1) 1 3 \ 

dp‘^ 

pdp /32 

\ 3(1 -b sin^ po) sin^ 2po ) 


+u{i3)-e^{L) tl;L0)=O, (7) 


and 


1 92 

W)d{¥y 


+ w{^^) - 


X(<^±) = 0, (8) 


where ^'J(/3,(/)) = N^ipL0)N^{x{^^)±x{^ ))/V 2 , with 
Np and being normalization factors, while {(3^) is the 

average of /32 over and = ^i^iL) + e^. 

On the above the following comments apply: 

a) (p = 0 corresponds to quadrupole deformation alone, 
while (p = ±t:I2 corresponds to octupole deformation 
alone. 

b) Because of the two steep oscillators involved, p re¬ 
mains close to 3zpQ and, therefore, the relative amount of 
quadrupole and octupole deformation remains constant. 

Some details of the derivation of the Hamiltonian, the 
method of solution, and the comparison of the present 
model to other approaches are given in Appendices Dl- 
D4. 


B. The d part of the spectrum 

Eq. 0 is exactly soluble [53l [54] in the case of the 
Davidson potentials [22] 




(9) 


in which the eigenfunctions are Laguerre polynomials 




with 


r(n^ + 0^ + 1) 




while the energy eigenvalues are given by 

Enfi,L = 2n^ -|- Oi -|- 1 = 2np + 1 

-b 


/ L(£-H) 3 

3(1 -b sin^ po) sin2 2pQ 


+ /9o"- (12) 


Eq. 0 is also exactly soluble in the case of an infinite 
well potential in /3, in which the eigenfunctions are Bessel 
functions. This solution has been worked out in Ref. m- 


C. The p part of the spectrum 

Eq. 0 with the potential corresponding to two har¬ 
monic oscillators centered at ± 0 o 

w{P^) = ^cipTpo)'^ = ^c{p^f, p^ = PtPo, (13) 
has been solved in Ref. [15]. The energy eigenvalues are 


^'l> = \l^ + ^ 710 = 0,1,2,... (14) 


m 


where is the number of quanta in the p degree of free¬ 
dom, while the eigenfunctions are Hermite polynomials 


H, 




b = 


c(/32) 


1/4 


, (15) 


with normalization constant Ab,, = . / • —r. 

?■ Y v’’'2 

The total energy in the present model is then 
E{np, L, po, n^) = Eq + CiE^^^l + C2n^, 




In what follows only bands with = 0 will be consid¬ 
ered. 

As an example of the dependence of the energy levels 
on the free parameters Po and po in Eq. (12), the energy 


3(1 -b sin2 Pq) sin2 2pQ 


levels with = 0 and L = 2 are shown in Fig. 1. 
Smooth variation with both parameters is seen. It is 
worth remarking that very slight dependence on pQ is 
observed between tt/S and Stt/S, in agreement with the 
findings of Ref. [15]. This observation (partly) justifies 
a posteriori the adiabatic approximation used in relation 
to the p degree of freedom, described in Appendix D3. 
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FIG. 1: (Color online) Excitation energies [Eq. (12l] 

for 71/3 = 0 and L = 2 are shown as a function of the tree 
parameters /3o and 4>o- All quantities shown are dimensionless. 
See subsection II.C for further discussion. 


D. B{EL) transition rates 


The electric quadrupole and octupole operators are 
(Eq. (6-63) of [J) 


= t2/32<i(0) = £/3cos<(>2?("i(0), 

rf 3) = ^ (17) 

with 


3Ze 3Ze ^3 

t2 = h = 

47r 47r 


(18) 


where R is the effective radius of the nucleus, while the 
electric dipole operator reads m 


rp{El) 




B 


\J B2B3 


52 sin 2()) (1) 

m./i 




(19) 


The total wave function reads 


$ 


L,M,q 


{P,cj7,e) = (/32/33)-'/"iV^F4(/3)iV^ 


± Xn^i4> )) 2L+1 


y /2 


327r2 


-(l±(-l)^)I?o^^^(d), 


( 20 ) 


where q stands for the rest of the quantum numbers (n^, 
n^). Since in what follows only bands with = 0 will 
be considered, in the remainder of the paper, as well as 
in the Appendices, we simplify the notation by using n 
instead of n,g. As a result, in what follows, rii and Uf 
indicate the initial and final values of np, while q/ and qf 
denote rii and Uf respectively, with = 0 at all times. 


B{EL) transition rates are given by 


B{EL\ L/qi —7 Ljqj) 


\{Ljqj\\T<^^R\\Uq,)\^ 

(2L. + 1) 


( 21 ) 


where the reduced matrix element is obtained through 
the Wigner-Eckart theorem 


( 22 ) 


In Eq. (21) the integration over the Euler angles 9 
involves standard integrals over three Wigner functions 
calculated in Appendix B, while the rest of the integra¬ 
tions are performed over / J where the /Sf, 

Z?! factors come from the volume element and cancel with 
the first factor of Eq. ( [^ . Using Eqs. (0 and ([^, 
as well as the relevant Jacobian, one finds ^p to con¬ 
stant factors) that the integration is over J f (i/32d/33 = 
-f ^dPdcj). 


Relevant integrals over (p are calculated in Appendix 
A, while integrals over jd are determined in Appendix 
C. The final results for matrix elements are summarized 
here. 

Matrix elements of between positive parity levels 
within the ground state band (rii = nf = 0) read 


{Lfqf\\T^'^\\Liqi)s^S — h 

y/2L, + l(Li2L/|000) 


B _ 1 cos </>o + e ^ 

-e 4b^ 




r(f+ ¥ + f) 


v/r(a, + l)r(ay + l)' 


(23) 


Matrix elements of from a positive parity level of 
the first excited band {rii = 1) to a positive parity level 
of the ground state band (n/ =0) are 


= t 2 \ 

V ^2 

_ 1 cos (in -|-, - 

® —Z 7572 —■\/2Ai -I- l(Li2I/y|000) 

r(f+ ¥ + i)r(f-^ + ^) 

v/r(a. + 2)r(ay + i)r - i) ’ 


(24) 


Matrix elements of between negative parity lev¬ 

els within the lowest band (n^ = n/ = 0) have the form 


IT Ilv.i52||r \ ^ _^COS(io -e ^ 

k^SdsWT =t2y^e 4 . 


^J2L, + l(Li2L/|000) 


r(^ + ¥ + f) 

v/r(a, + l)r(a/ + l)' 


(25) 


Matrix elements of between a positive parity 

level of the ground state band and a negative parity level 
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of the lowest band, or vice versa, read 

__ ^ 

B e 4i>2 sin^o 


{Lfqf\\T^^\\L,q,)=t,^ 


B 3 


v/2l7Ti 


r f^ _i_ it _i_ S') 

(L, 3 Lf| 000 ) , ^ . ( 26 ) 

^ Vr(«* + i)r(a/ + i) ^ ^ 

Matrix elements of between a positive parity 

level of the ground state band and a negative parity level 
of the lowest band, or vice versa, are 


or vice versa, read 


_ 1 • 2 / 

, , o B e ^ sin ©n 

B{E3;L,q,^ Lfqf)=tl- 


’B 3 (1 - 

fp I i/ 3')]2 

(Li3Lf|000)2-LA2—- dL, 

^ ^ r(a, + l)r(a/ + l) 


(31) 


B{El)s between a positive parity level of the ground 
state band and a negative parity level of the lowest band. 


or vice versa, are 


1 


{Lfqf\\T^ \\Liqi)s^S = 


B 


e ^ sin 2 (j)o 


VB 2 B 3 


p (^ _l_ it _l_ 2 ) 

a/ 2L, + l(LaLHOOO)^ ^ ^ ^ =■ (27) 

^ Vr(«* + i)r(a/ + i) ^ ^ 

The final results for B{EL)s are summarized here. 
B{E2)s between positive parity levels within the 
ground state band {rii = Uf = 0) read 


B{E 2 -Liq^ Lfqf )s^s = ^2-5“ 

1 {coscpo + e“^^‘^o)2 
e 2 i ,2 


(1 + 

{Li2Lf\000)‘ 


„-h2 02.2 


[r(^ 


i/ 


+ §)]' 


r(ai + i)r(a/ + 1 ) 


• (28) 


B{E2)s from a positive parity level of the first excited 
band (rii = 1) to a positive parity level of the ground 
state band (n/ = 0) are 


B{E2; Liq, Lfqf)s^s = ^2“^ 

B 2 

= -b"0o'|2 


(cos 00 


[r(^ + ^ + i)r(^-f+ 1)]^ 

r(a, + 2)r(a/ + i) [rif-'^-DY 


(29) 


B{E2)s between negative parity levels within the low¬ 
est band (n^ = nj =0) have the form 


B(^E 2 ,^ LiQi >■ = ^2^“ 

B 2 

_ 1 (cos (j)o 


2 \ r' 


(1 


— e ^^“^ 0)2 


[p i/ _1_ 3')12 

B{E3)s between a positive parity level of the ground 
state band and a negative parity level of the lowest band. 


^ , 1 2 73^ e b^sin^20o 

B{El-L,q, ^ Lfqf) = 

(^*1^/1000) r(a, + i)r(ay + i)- 

In Ref. [15] it has been pointed out that B{El)s for 
transitions from positive parity levels L of the ground 
state band to negative parity levels L — 1 of the lowest 
band and B{El)s for transitions in the opposite direc¬ 
tion, i.e. from negative parity levels L — 1 of the lowest 
band to positive parity levels L of the ground state band, 
should be of the same order, as seen experimentally. This 
condition is clearly fulfilled by Eq. (32). 


III. NUMERICAL RESULTS 
A. Spectra 

In the spectra of the AQOA model with Davidson po¬ 
tential (AQOA-D) only the parameters 0o and (3o play 
an essential role, as seen in Eq. (121, while the quantities 
Eq, Cl, and C 2 of Eq. (16) do not enter, if we consider 


only bands with n,p = 0 and normalize all energies to 
that of the Hrst excited state, E{2Y). The parameters of 
rms fits, using the quality measure 


a = 


l Etimexp) - E^hW 

(n-l)E(2+)2 


(33) 


to the spectra of the Ra and Th isotopes lying at the 
border between the regions of octupole deformation and 
octupole vibrations, as well as within the former region, 
are shown in Table I, while in Table II the relevant spec¬ 
tra are shown. For ^^®Ra and ^^®Th, the predictions of 
the original one-parameter (0o) AQOA with an infinite 
well potential (AQOA-IW), applicable at the border be¬ 
tween the regions of octupole deformation and octupole 
vibrations, are shown for comparison. The following com¬ 
ments apply. 

a) Good agreement between the theoretical predictions 
of AQOA-D and experimental data is obtained up to high 
angular momenta, both in the ground state band and in 
the negative parity band. The predictions for the 1“ and 
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TABLE I: Parameters 4>o and /So of the AQOA model with Davidson potential (AQOA-D) obtained from rms fits to experimental 
spectra of ^^^Ra [2l[8], ^^"^Ra [ZlE], ^^®Ra [ZllS], ^^"^Th [24] . and 2^®Th The experimental -R4/2 = j5(4^)/i7(2^) ratios are 
also shown. The angular momenta of the highest levels of the ground state, /3 and negative parity bands included in the rms fit 
are labelled by Lg, Lp, and Lo respectively, while n indicates the total number of levels involved in the fit and a is the quality 
measure of Eq. ( |33[ ). The theoretical predictions are obtained from the formulae mentioned in subsec. II.B . See subsec. III.A 
for further discussion. For completeness, the parameter b, and the parameter ratios t 3 ,/t 2 , and B 2 /B 2 ,, appearing in the fitting 
of electric transition rates, are shown, wherever known. All quantities shown are dimensionless, exept t3/t2, which is given in 
fm. See subsec. III.B for further discussion. 


nucleus A4/2 

00 

00 

L, 

Lo 

n 

(J 

b 

0/0 (fm) B 2 IB 3 

^^""Ra 2.715 

62.7 

0.00 

20 

0 

19 

20 

0.917 


7.266 


^^^Ra 2.970 

41.9 

1.80 

28 

0 

27 

28 

1.351 

1.836 

7.288 

0.65 

22®Ra 3.127 

24.7 

2.14 

28 

0 

27 

28 

1.360 


7.309 


^^^Th 2.896 

67.9 

0.79 

18 


17 

17 

0.843 


7.288 


2^®Th 3.136 

9.5 

0.94 

20 

0 

19 

20 

0.994 


7.309 



TABLE II: Comparison of theoretical predictions of the AQOA model with Davidson potential (AQOA-D, columns labelled 
as D) and with an infinite well potential (AQOA-IW, columns labelled as IW) to experimental data [normalized to i?(2(^)] 
of ^^^Ra [71 H], ^^"^Ra [71[H], ^^®Ra [71 |H], ^^'*Th |^, and AQOA-D parameters are shown in Table I, while 

the AQOA-IW predictions have been taken from Ref. |15| . The theoretical predictions for AQOA-D are obtained from the 
formulae mentioned in subsec. II.B . See subsec. III.A for further discussion. 







''''®Ra 

''''®Ra 

''''®Ra 

^""^Th 

22 irj.^ 




L'^ 

exp. 

D 

exp. 

D 

exp. 

D 

IW 

exp. 

D 

exp. 

D 

IW 

4+ 

2.72 

3.00 

2.97 

3.17 

3.13 

3.22 

3.09 

2.90 

3.09 

3.14 

3.22 

3.12 

6+ 

4.95 

5.59 

5.68 

6.21 

6.16 

6.45 

5.99 

5.45 

5.90 

6.20 

6.44 

6.10 

8+ 

7.58 

8.49 

8.94 

9.87 

9.89 

10.45 

9.56 

8.50 

9.17 

10.00 

10.42 

9.78 

10+ 

10.55 

11.58 

12.66 

13.94 

14.19 

15.02 

13.71 

11.97 

12.71 

14.41 

14.97 

14.08 

12+ 

13.82 

14.77 

16.74 

18.30 

18.93 

20.00 

18.42 

15.80 

16.42 

19.32 

19.93 

18.96 

14+ 

17.39 

18.04 

21.17 

22.85 

24.06 

25.31 

23.64 

19.97 

20.25 

24.68 

25.20 

24.38 

16+ 

21.21 

21.36 

25.90 

27.55 

29.52 

30.84 

29.38 

24.44 

24.16 

30.41 

30.69 

30.34 

18+ 

25.28 

24.70 

30.92 

32.34 

35.30 

36.54 

35.61 

29.20 

28.13 

36.50 

36.35 

36.81 

20+ 

29.57 

28.07 

36.22 

37.22 

41.38 

42.38 

42.33 



42.90 

42.14 

43.80 

22+ 



41.74 

42.15 

47.75 

48.32 

49.54 






24+ 



47.48 

47.13 

54.44 

54.34 

57.22 






26+ 



53.41 

52.14 

61.42 

60.43 

65.38 






28+ 



59.54 

57.18 

68.70 

66.57 

74.01 






0+ 

8.23 

8.06 

10.86 

10.90 

12.19 

12.21 

11.23 


9.91 

11.18 

11.31 

12.41 

1“ 

2.18 

0.35 

2.56 

0.34 

3.75 

0.34 

0.34 

2.56 

0.34 

3.19 

0.34 

0.34 

3- 

2.85 

1.90 

3.44 

1.95 

4.75 

1.97 

1.93 

3.11 

1.93 

4.26 

1.97 

1.94 

5" 

4.26 

4.24 

5.13 

4.60 

6.60 

4.72 

4.45 

4.74 

4.43 

6.24 

4.72 

4.51 

7" 

6.33 

7.01 

7.59 

7.98 

9.26 

8.36 

7.70 

7.13 

7.49 

9.11 

8.35 

7.86 

9" 

8.92 

10.02 

10.73 

11.86 

12.68 

12.67 

11.57 

10.17 

10.91 

12.79 

12.64 

11.86 

11“ 

11.97 

13.17 

14.46 

16.09 

16.74 

17.47 

16.00 

13.73 

14.55 

17.15 

17.41 

16.45 

13" 

15.38 

16.40 

18.68 

20.56 

21.39 

22.63 

20.96 

17.72 

18.32 

22.11 

22.53 

21.61 

15" 

19.11 

19.69 

23.31 

25.18 

26.54 

28.05 

26.45 

22.07 

22.20 

27.55 

27.92 

27.30 

17- 

23.11 

23.03 

28.27 

29.93 

32.13 

33.67 

32.43 

26.71 

26.14 

33.42 

33.50 

33.51 

19” 

27.35 

26.39 

33.51 

34.77 

38.10 

39.44 

38.91 



39.63 

39.23 

40.24 

21" 



38.99 

39.68 

44.41 

45.34 

45.87 






23" 



44.67 

44.63 

51.03 

51.32 

53.32 






25" 



50.55 

49.63 

57.93 

57.38 

61.24 






27" 



56.60 

54.66 

65.08 

63.49 

69.64 







3 states are poor, since no finite barrier is used in the 
phi potential [5M57] 

b) When moving from the border region towards the 
interior of the region of octupole deformation, the pa¬ 
rameter 00 increases, in agreement with an increasing 
role of the octupole deformation (which is proportional 


to sin0o)j while the parameter /3o decreases. In paral¬ 
lel, a decreasing role of the quadrupole deformation is 
revealed by the decreasing Ri /2 ratios. 

c) The agreement between the predictions of AQOA- 
IW and the data is comparable to that of AQOA-D, ex¬ 
cept at high angular momenta, where AQOA-D is closer 
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TABLE III: Matrix elements of electric transitions in 
The experimental data, in units of e fm, e fm^, e fm^ for El, 
E2, E3 respectively, have been taken from Ref. |21| . while the 
theoretical predictions have been obtained using the formulae 
of subsec. II. D in the way described in subsec. III.B, where 
further discussion is given. 


m.e 



exp. 

th. 

70^ 

\E2\ 

|2+) 

199±3 

196 

(2+ 

\E2\ 

|4+) 

315±6 

323 

(4+ 

\E2\ 

|6+) 

405±15 

426 

(6+ 

\E2\ 

|8+) 

500±60 

525 

(1- 

\E2 

|3-) 

230±11 

236 

(3- 

\E2 

|5-) 

410±60 

334 

(0+ 

\E2\ 

|2+) 

23±4 

36 

(0+ 

\E3\ 

|3-) 

940±30 

1006 

(2+ 

E3 

|1-) 

1370±140 

1137 

(2+ 

\E3\ 

|3-) 

<4000 

1176 

(2+ 

E3 

|5-) 

1410±190 

1594 

(0+ 

El 

|1-) 

<0.018 

0.013 

(2+ 

El 

|1-) 

<0.03 

0.018 

(2+ 

El 

|3-) 

0.026±0.005 0.023 

(4+ 

El 

|5-) 

0.030±0.010 0.032 

(6+ 

El| 

|7-) 

<0.10 

0.042 


to the data, due to the term /?q appearing in Eq. 
which moderates the increase of the energy with L. 


( 121 , 


B. Transitions 


With the parameters Po determined from the 

spectra, we now turn attention to electromagnetic tran¬ 
sition rates, following the procedure described below. 


a) From Eqs. (28) and (29) it is clear that ratios of 
B{E2)s^sS involve only the parameters 4 >o and /3o, thus 
they are already fixed. The same holds separately for 
ratios of B{E2)a^as, or B{E3)s, or B{El)s, as seen 
from Eqs. ([^, (31), and ([^ respectively. 


b) Ratios of B{E2)s^sS over B{E2)a^aS involve in 
addition the parameter &, which can then be determined 
from such ratios, as seen from Eqs. (p^ and (30). 


c) Ratios of B{E3)s over B{E2)s involve in addition 
the ratio t 3 jt 2 , which can be determined from the nuclear 
radius (see subsec. III.B.2), and the ratio R 2 /S 3 , which 
can be determined from the B{E3)/B{E2) ratios, as seen 
from Eqs. (31) and (28). 


This procedure can be tested against the recently mea¬ 
sured transition matrix elements of ^^^Ra |21j , shown in 
Table III. 


1. E2 transitions 


The ratio of any E2a^a matrix element [Eq. (25)] 
over any E2s^s matrix element [Eq. (23)] contains the 


ratio 

^ ^ (cos0o-a)(l-ba) ^ ^ 
{cos(j)o -I- a)(l — o) ’ 


(34) 


In the case of ^^^Ra, the ratios of experimental matrix 

elements |+^, 4 + 1 ( 2 + . 6 +^ 4 + : 1 +^1+ lead to r = 
0.852, 0.887, 0.909, 0.908, i.e., to an average value of 
0.889 . 

Solving for a one obtains the quadratic equation 
(r — l)a^ -I- (r-l- l)(cos(()o — l)a-|- (1 — r) cos(j)o = 0, (35) 


having the solution 


[-(r -b l)(cos^o - 1) 


2(r-l) 

±-\/(rT l)2(cos())o — 1)^ -|- 4(r — 1)2 cos4>o 


Since from Eq. (34) one has 


52 = - 


In a 




(36) 


(37) 


In a has to be negative for b to be real. Since r < 1, we 
see that in Eq. (36) only the negative sign is allowed 
in order to have a > 0, leading in the case of ^^“^Ra to 
a = 0.1647 and b = 1.836 . 

One can then keep 


E = t2 \J B/B 2 


(38) 


as an overall constant for all E2 transition matrix ele¬ 
ments, and determine it through rms fitting to the ex¬ 
perimental data of the transitions with = n/ = 0, 
obtaining F = 127.20 and the E2 predictions reported 
in Table III. Then from Eq. (24) one can calculate also 


the transitions with = 1, n/ = 0, one of which is also 
reported in Table III. 


2. E3 transitions 


For the coefficients t 2 and tg one can use Eq. (18), 
leading to 


^3 _ n 

where R is the nuclear radius, given by 

R = rg = 1.2 fm. 


(39) 


(40) 


with A being the mass number of the nucleus. Thus in 
the case of one has R — 7.2878 fm. 

One can then determine the ratio ^ from any ratio of 
E3 matrix element over E2 matrix element, since each 
of these ratios contains the quantity { 1 ^/ 12 )^/B 2 /B-^, as 
seen from Eqs. (26), (23), (25). In the case of ^^^Ra, 


six E2 matrix elements and three E3 matrix elements 


































are known. Considering all 18 pos sible ratios, we get 
an average value of {H/t 2 )\/B 2 /= 5.876, leading to 
B 2 /B 3 = 0.65. 

By now the ES transition matrix elements have been 
completely determined. The overall constant 

F' = H^BIB^ (41) 

appearing in this case is connected to the overall constant 
F through 


F' _ ts I B2 _ I B2 

F ^2 V ^3 V ^3 ’ 


(42) 


leading to F' = 747.47 and to the E3 matrix elements 
given in Table III. 


3. El transitions 


In the case of El matrix elements, the quantity 


1 B e ^ sin 2^0 

2 ^ ^B2B3 y/l- 


(43) 


can be treated as an overall constant, determined in the 
case of ^^^Ra by rms fitting to the two known transitions 
to be F” = 2.676 10“^ and providing the predictions 
given in Table III. 


IV. CONCLUSIONS 

The analytic quadrupole octupole axially symmetric 
model with an infinite well potential (AQOA-IW) had 
successfully predicted the border between the regions 
of octupole deformation and octupole vibrations in the 
light actinides, identifying ^^®Ra and ^^®Th as border 
nuclei with heavier isotopes corresponding to oc¬ 
tupole vibrations and lighter isotopes exhibiting octupole 
deformation. The AQOA-IW model involved only one 
free parameter, (j)Q, expressing the relative presence of 
quadrupole vs. octupole deformation, while a parameter- 
free version has also been developed later [59]. 

In the present work, the infinite well potential is sub¬ 
stituted by a Davidson potential, resulting in the AQOA- 
D model, which is able to deviate from the border line 
into the region of octupole deformation. This is achieved 
through the extra parameter /3o, the position of the min¬ 
imum of the Davidson potential, which is increasing with 
increasing R 4/2 ratios, as it is known from its use in the 
description of quadrupole deformed nuclei [53] . 

Within the AQOA-D model, analytic expressions for 
energy spectra and B(EI), B(E2), B(E3) transition rates 
are derived. Then the following path is taken. 

a) The spectra of and 224,226rpj^ [normalized 

to E[ 2 f )] are well reproduced in terms of the above men¬ 
tioned two parameters ipo and Pq. 


b) The parameter b, related to the harmonic oscillator 
potential used in the (p degree of freedom, can be deter¬ 
mined from the ratio of any E 2 matrix element between 
negative parity states over any E 2 matrix element be¬ 
tween positive parity states, fixing the determination of 
all E 2 transitions up to an overall scale factor. 

c) The ratio of mass parameters i?2/7?3 can be deter¬ 
mined from the ratio of any E3 matrix element over any 
E 2 matrix element, while the ratio of transition coeffi¬ 
cients t 2 /t 3 is fixed by the nuclear radius. As a result, 
the determination of all E 2 and E3 transitions is fixed, 
without any additional overall scale factor. 

d) El transitions are also fixed, up to another scale 
factor. 

The recently measured B{EL) transition rates of ^^^Ra 
[21] . presenting stable octupole deformation, provide a 
successful test for the model. It is clear that for other 
nuclei, the minimum set of data needed includes 

a) A few energy levels of both positive and negative 
parity, from which the parameters (po and Pq can be de¬ 
termined. 

b) At least one E 2 transition between positive par¬ 
ity states and one E 2 transition between negative parity 
states, from which the parameter b can be determined. 

c) At least one ES transition, from which, in combina¬ 
tion with the E 2 transitions of b), the parameter ratio 
B 2 /B 3 can be determined. 

From these pieces of data 

a) The spectrum (leaving out the 7 bands) is deter¬ 
mined up to an overall scale factor. 

b) All relevant E2 and E3 transitions are determined 
up to an overall scale factor. 

c) All relevant El transitions are determined up to 
another overall scale factor. 

It is of interest to apply the present model in the ac¬ 
tinides close to ^^°Pu, in which a second order shape 
phase transition from octupole-nondeformed to octupole- 
deformed shapes has been recently found ED], while oc¬ 
tupole bands have been described [6T] using supersym¬ 
metric quantum mechanics. The light rare earths, in 
which octupole bands have been considered recently both 
by the Bizzeti and Bizzeti-Sona approach [62] and within 
density functional theory [63], are also of special interest. 
A successful application of the AQOA-IW model to ^^®Nd 
has already been given in Ref. [64] . 
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Appendix A. tf) integrals 


For antisymmetric states one has 


Since we confine ourselves to states with = 0, this 
quantum number is omitted in the notation of the wave 
functions, which then carry only the subscript i (/) for 
the initial (final) state. 

For symmetric states one has 


Xsi^)=Ns 


72 


while for antisymmetric states one has 

x{^'^)-x{^~) 


XaW=Na^ 


72 


(44) 


(45) 


where Ns and Na are normalization factors and, accord¬ 
ing to Eq. ([T^, 


x(</'+) = 



) = 7 ^ 





AO. Normalization 


(46) 


^ = / XAXAdck 


> —oo 

1 r 


1 

: J — OO 

pOO 

- / xi(l>^)xi(t^~)d(l) 


1 b 

27^' 


o-bNi 


' —oo 
/*oo 


^_b2 02_2t,2 


20f 


“ ~1=^ 
\/tt 


b -b^^l 


e-^ '^^dcj), (50) 


leading in the same way as above to 
1 


Nl 


= l-e-^^o. 


(51) 


Al. B{E2)s 

The transition operator for B{E2)s contains /32 = 
/3 cos (j). 

For B{E2)s between symmetric states one has 


For symmetric states one has 


^ = / XsXsdcj) 


' —OO 

a 


1 r°° 

{xi(t>'^)fd<j)+- {x{(l)~)fd(j) 

: ^ J — OO 


xi(l>~^)xi^ )d(t) 


1 b 


' —OO 

u2,2 


2v^ 

1 b 

27^' 


^-b-‘<p^+2b- 




' —OO 
pOO 


-26- 


+ r e-^"*"d(j). (47) 

V J — OO 


Using Eq. (|84|) of Appendix A4 we see that the integrals 

(48) 


appearing here are of the form 

pOO 

/ -b^4>^±2h 




leading to 


1 1,2 J 2 

1172 = 1 + e" ■ 


t(.E2) 


/ Xs cos (j)Xsd4> 

-OO 

/y-2 poo 

= ^y {x{(l}^)f cos (l)d(l) 

]Sr2 fOO 

+ ^J (x(^”))^cos0# 

/ OO 

x(<(>’^)x(«<>")cos (t>d(t) 

-OO 


A 

2 7^ 


2 7^ 


-b^(p^-\-2b^4>o(p 


-b^ 4>^-2b^ M 


COS (pdcj) 


cos (pdcj) 


9 b 

+ '^0 / 

7^ j-i 


s ^ cos 4>d<j). (52) 


Using Eq. (81) of Appendix A4 we see that the integrals 
appearing here are of the form 

[ e~^ ^ cos (j>d(j) = ‘l>o-Ab^ cos(j)o, (53) 

J-oo b 

leading to 


2^2 

-^0 


r{E2) _1.^ COS ((iQ + e 

A4s = e ^ 


1 + 


(49) 


(54) 
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For B{E2)s between antisymmetric states one has 


For B{E3)s between symmetric states one has 


AE2) 

'(p,A^A 


Xa cos (fiXAdtj) 

) 

Nl 


= ^ / cos (j)d(j) 

^ J — OO 

foo 

+ ^y cos (j)d(t) 

/ OO 

x(^^)x(</>”)cos (l)d(j) 

-OO 

-b"0^-2b^,^o<^(,OS^# 




— OO 
2^2 


2 0r 

2 i/tt 




_iV2_^e-''^o 


e ^ cos 4)d(f), (55) 


leading in the same way as above to 


AE 2 ) _ 1 COS(/io -e ‘ 

h,A^A = C 


(56) 


For B{E2)s between symmetric and antisymmetric 
states one has 


AE2) 

^4>,S^A 


Xs COS (j)XAd(l) 

30 

^sNa rt ,,+,,2 ,,, 

(x(</> )) cos0d(/) 

(x(<(>"))^cos 4>d(j) 


2 

NsNa 


' —OO 

/*oo 


NsNa b _j,2 ,2 


e / g-6"0=+2h"0o<^ 


2 0r 

2 


leading to 


r(-E2) _ Q 

^d,,S^A — ^• 


In the same way one also finds 


r(-E2) _ Q 


COS (j)d4> 

cos (j)d(j), (57) 

(58) 

(59) 


A2. B(F3)s 

The transition operator for B{E3)s contains /Sa 
P sinc^. 


AE3) 

>,S^S 


/ Xs sin (j)Xsd(j> 

J —OO 

N "2 /-OO 

= ^y (x(</>’^))^sin(/)d(;i 

N ‘2 fOO 

+ ^J smpdp 

/ OO 

x(</>'^)x(<('")sin (t>d(j) 

-OO 


2 


' —C30 
rOQ 


0-62 02+262 00 ,^ 


-62 02-262 000 


sin (j)d(j) 


sin (pdcj) 


2 i/tt 

+'^0 [ e~^sin (j>d(j). (60) 

V ^ J —OO 


Using Eq. (83) of Appendix A4 we see that the integrals 
appearing here are of the form 

J-oo b 

leading to 


E j, 2^2 _ 1 _ 

4|>2 sm ( 


(61) 


r(S3) 


= 0 . 


p,s^s — (®2) 

For B{E3)s between antisymmetric states one has 


AE3) 

'<I>,A^A 


/ XA sin (f)XAd(f) 

J —OO 

t\j2 fOO 

= ^ / {x{<l}~^)f sinpdcl) 

^ J — OO 
/V^ 1*00 

+ ^ / (x(<(’"))^sin<()(i(() 

J —OO 

/ OO 

x(</>’^)x(<(’")sin pdp 

-OO 


^jLg-b'^'Pl [°° g-6202+26= 000 
2 J-oo 

Xo b _iAaA f°° -6=0=-2b=0o0 


bS UUg-6 


sin (j)d(j) 


sin pdcj) 


2 i/tt 

E ^OO 

— Ns—^e~^‘^° / 6“^*^ sin (j)d4>, (63) 

Vd-oo 


leading in the same way as above to 


t{E3) _ Q 


(64) 


For B(E3)s between symmetric and antisymmetric 
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states one has 


/ Xssin4>XAdc^ 

^ —oo 


NsNa 

2 

NsNa 


— OO 
POO 


sin 

(x(<^"))^sin (j)d<j3 


2 ^y^: 

NsNa b _j,2^2 

2 ^/7r 


' —OO 
^OO 


/ oo 

g-b 0 -2b ^05^ 

-oo 


leading to 


/ 


(BS) _ e ^ sin (^q 


( 66 ) 


= (67) 


In the same way one finds 

[£3) 

t>,A^S ~ ^<j>,S^A- 


A3. B(£;i)s 

The transition operator for B{El)s contains P 2 P 3 = 
p cos PP sin P = P"^ Sin 2 (ji ^ 

For B{El)s between symmetric states one has 

1 

2 J_^Xssin2pXsdp 

fa2 foo 

= ^ / ixiP'^))'^ Sin2pdp 

^ J —OO 
N2. fOO 

+ ^ ixiP~)f sin2pdp 

^ J — OO 
N"2 r°° 

+ ^ x{P~^)x{P~)sin2pdp 


' —OO 

^00 

=e ■' / e 

> —00 




M jLp-b"0o -b"<^.= -2b^0o0 

4 7-00 


.2^2 I 

2 


sin 2pdp 
sin 2pdp 

00 

e“^ ^ sin20£i(/). (68) 


Using Eq. (83) of Appendix A4 we see that the integrals 


appearing here are of the form 


[ e ^ sin2pdp = <>2 sin2(()o, 

7-00 b 

(69) 

leading to 


j(si) _o 


For _B(Ul)s between antisymmetric states one has 

POO 


I 


(El) 

},A^A 


1 


Xa sin 2pXAdp 


+ 


r°° 

/ ixiP'^yf sm2pdp 

^ J — 00 


(x(^ ))^ sin 20^0 


/ —00 
/•oo 


TV? Z*'^ 

-^ / X(0^)x(7' )sm2(^d(() 

^ t/ —00 

= r e-'>"^"+^^"^-*sin2pdp 

4 7-00 


4 

_ 

2 




sin 2pdp 


A e-b"<?ig 


e ^ sin20£i^, (71) 


leading in the same way as above to 


t(E1) 

t.A^A 


(72) 

For B(El)s between symmetric and antisymmetric 
states one has 


r(Bl) 

S^A 


1 


Xs sin 2pXAdp 


NsNa 

4 

NsNa 


' —00 

pOO 


{x{p'^)f Sin2pdp 
ix{P~)f sin2pdp 


_ NsNa b 52^2 

4 

NsNa b ^ 2.2 


„-b'^<P+2b^<t>o<t> 


sin 2pdp 


4 Vtt 
leading to 


/ OO 

e-''"'^2-262000 (73) 

-00 


AEi) _ e sin2(/)o 

~ 2 V 1 - e-2'>% 

In the same way one finds 


/ 


.(£1) ^ jiEl) 


>,A^S 


',S^A- 


(74) 


(75) 


A4. Useful integrals 
We know that (Eq. 3.897.2 of Ref. [55] 1 

POO 

/ cos bxdx = 

Jo 

I /tt 


4 V 6 


,(7-*b)U(4/3) 


|g(7+ib)V(4/3) 


2VP J 

1-$(1^ 


, (76) 


(70) 
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where Ref3 >0, b > 0, and $( 0 ;) is the error function, 
having the property 

$(-x) = -$(x). (77) 


Changing the variable into y = —x, Eq. (76 1 takes the 
form 


p — OO 

/ cos bydy = 

Jo 


1 /tt 


dV 6 


(7-*6)V(4/3) 

(7+*t.)V(4/3) 


V 2v^ J\ 

7 + 


n-P> J., (78) 

Changing the symbol y into x and letting 7 —)■ — 7 , one 
then gets 

/-o 


1 /TT 


dV 6 ^ ® 


cos bxdx = 


(-7-*f>)V(4/3) 


1 - $ 


+e 


(-7+ib)V(4/3) 


1 - $ 


—7 — ib 
2v^ 
—7 + z 6 

2V? 


■ (79) 


Taking into account the property ( ff^ , Eq. ( [ 7 ^ takes 
the form 


cos bxdx = 


1 /TT 


dV 6 


+e 


(7+2f>)V(4/3) 

(7-2b)V(4/3) 




l + d> 


7 — i 6 \ 


(80) 


Adding Eqs. (76) and (80), we get 


/ e ^ ^ cos bxdx —I e + e 

^ — CXD \ P \ 

[tt 76 

= V^' " <“> 

In a similar way, starting from the integral (Eq. 3.897.1 
of Ref. [55] ! 




_ 1<j gh-ibf/m 


—e 


{7+2f>)V(4^) 


- Kw ) 


(82) 


where i?e/3 >0, & > 0, we get 


e smbxdx =- le - e 

1 2 y p \ 

pK y^-b^ . 76 

= “7' " 


For normalization purposes the integral 


J ^-(ax^+bx+c)^^ ^ y^g(b=-4ac)/(4a) 


suffices. 


Appendix B. 0 integrals 

Integrals over 9 involve three Wigner functions and can 
be calculated using Eq. (d.6.2) of Ref. [ 66 ] 

/ (0)17^1(0) 

= 8^" { f f f 1 • (85) 

\ ki k 2 ks J \mi m 2 m 3 J 

Using the relation between 3-j symbols and Clebsch Gor- 
dan coefficients (3.7.3) of [ 66 ] 

ji j2 j3 \ ^ ^-^y!:!::: i(^j^j^j^irn,m2-m3), 
mi m 2 m 3 J 72 ja + 1 ''4442:431 2 3 J, 

( 86 ) 

and the relation for conjugate Wigner functions (d.2.7) 
of [ 66 ] 

= (-l)'=-'”27i'L™(^), (87) 

one obtains 

/ milt. m-m 

8^2 

(jlj 2 j 3 |fclfc 2 - fc3)(jlj2j3|TOlTO2 - TO3), (88) 


2J3 + 1 

which by replacing m3 (fca) by —m3 {—k 3 ) can be rewrit¬ 
ten as 

8^2 

= tt- ——r(U-?2j3|fciA2fc3)(jij2j3|™im2m3). (89) 

2J3 + 1 

Bl. B{E2)s 

In this case the integral reads 

7^^ = / d^o^^-^[i ± im'PtM.mZid) 
^^0^[l±(-l)^/PoX(^)- (90) 


Using Eq. (89) this gives 

\L 


AE2) _ (1±(-1)^-)(1±(-1)^4) / 2L. + 1 

« d V 2 T/ + 1 

(L,2L/|000)(L,2Ly|M,pM/). (91) 
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From the (j) integrals we know that non-vanishing re¬ 
sults are obtained only in the 5” —>■ S' and A cases. 


In the S ^ S case the two factors in the rhs of Eq. (91) 


have the positive signs in the place of the double signs, 
thus allowing only even values of Li and Lf , resulting in 
a factor of 4 in the numerator. 

In the A —> A case the two factors in the rhs of Eq. (91) 


have the negative signs in the place of the double signs, 
thus allowing only odd values of Li and L/, resulting 
again in a factor of 4 in the numerator. 

As a consequence, in all cases the final reasult reads 


j{E2) _ / -I- 1 

« " V 2^/ + 1 


{L,2Lf\mA){Li2Lf\Mi^jLMf). (92) 


B2. B{E3)s 


The calculation parallels the one of the previous sub¬ 
section, the only difference being that the middle term, 

(3) 

coming from the transition operator, is I?g The result 
reads 


jiE3) ^ (1±(-1)^-)(1±(-1)^0 / 2T, + 1 
® 4 \j2Lf + l 

{Li3Lf\000){LiSLf\MinMf). (93) 


From the (j) integrals we know that non-vanishing re¬ 
sults are obtained only in the 5” —>■ A and A —>■ S' cases. 

In the S —)■ A case, Li is even and Lf is odd. The 
first factor in the rhs of Eq. (93) has the positive sign in 
the place of the double sign in front of the (—1)^’ term 
and the negative sign in the place of the double sign in 
front of the (—1)'^-'' term, resulting in a factor of 4 in the 
numerator. The same factor of 4 is obtained also in the 
A —>■ S case. Therefore in all cases the final result reads 


AE3) _ 

In - 


'2L, + 1 
2Lf + 1 


(L,3L/|000)(L,3L/|M,/xM/). (94) 


B3. B{El)s 


The calculation parallels the one of the previous sub¬ 
section, the only difference being that the middle term, 
coming from the transition operator, is The result 

reads 


-1) _ (1±(-1)^-)(1±(-1)^/) / 2L. + 1 

® 4 \j2Lf + l 

iLALf\000){LALf\M,AMf). (95) 


Erom the 4> integrals we know that non-vanishing re¬ 
sults are obtained only in the S' —>■ A and A —>• S cases. 

In the S —)■ A case, Li is even and Lf is odd. The 
first factor in the rhs of Eq. (95) has the positive sign in 


the place of the double sign in front of the (—1)^’ term 
and the negative sign in the place of the double sign in 
front of the (—1)-^^ term, resulting in a factor of 4 in the 
numerator. The same factor of 4 is obtained also in the 
A —7> S case. Therefore in all cases the final result reads 


4""'^ = J^^^{LdLf\000)iLdLf\M,AMf). (96) 


B4. Normalization 


Normalization of the 9 wave functions is guaranteed by 
the integral of Eq. (4.6.1) of Ref. [Ml 

(97) 

The normalization integral for any state reads 


= / d'^9 


2L + 1 
327r2 


[1±(-1)Wm(0) 


(98) 


Using Eq. (971 this gives 

^ H-(-1)2^±2(-1)^ 


= 1 , 


(99) 


since for symmetric states the positive sign appears in 
the place of the double sign and L is even, while for an¬ 
tisymmetric states the negative sign appears in the place 
of the double sign and L is odd. 


Appendix C. /3 integrals 

Cl. B{E2)s 

The transition operator contains a /3 factor, thus the 
integrals appearing in this case read 


l{E2) ^ 


with 


I 2(n,)! 2(n^)! ^^^ 2 )., ^ , 

r{ni + Qi + 1) r{nf + af + 1) 

( 100 ) 


nf) = / {P^)LZ iP^)dp. 

Jo 

(101) 

Using the substitution = x with dx = 2l3dl3, the inte¬ 
gral is written as 

1 POO ^ ^ 

I^^'^\ni,nf) = - / e“''x^ + ~ + 5L“((a:)L5(^(a;)dx. 

^ JO 

( 102 ) 
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Analytic expressions can be found for these integrals 
in the case in which one of the quantum numbers rii, Uf 
is zero. (For the case in which both quantum numbers 
m and Uf are non-zero, see Eq. (B5) of Ref. [ST].) We 
consider n/ = 0, since both the ground state band and 
the octupole band are characterized by this value. Then 
one has L“^^g(a:) = 1 and the integral is simplified into 


C2. B{E3)s 

The transition operator again contains a /? factor, thus 
the integrals appearing in this case are exactly the same 
as in the previous subsection 




e ^ ^ L^\[x)dx. (103) 


Integrals of this form are known to have the following 
analytic solution ([68], p. 463, Eq. (5)) 


^ j(e2 ) 


C3. B{El)s 


( 112 ) 


e ^^L^{cx)dx = {a), 

where (a)n is the Pochhammer symbol 

r(a -I- n) 


(a)n = 


r(a) 


(104) 


(105) 


The transition operator contains a factor, thus the 
integrals appearing in this case read 


j(El) ^ 


with 


2(n,)! 


T{ni -I- ttj -I- 1) r(n/ + af -\-l) 




( 113 ) 


“f 

2 


+ 75, i.e. 


By replacing c = 1, A = a^, a — 1 = ^ 
a = ^ ^ and applying the definition (1051 the 

result is 

l(E2)r QT) _ 1 r (Uj -t f - ^ - A) /3 

(106) 

In the simplest case of a transition between states with 
Hi = 0 and nj = 0, which will be eventually of major 
interest in the present work, one obviously has 




Qj a f 3 
2 2 2 


(107) 


Jo 

. (114) 

Using again the substitution = x with dx = 2/3d/3, the 
integral is written as 

1 a 

I^^^\ni,nf) = - / e““a:^ + ^+^L“'.(x)L5J^(x)dx. 

(115) 

For Uf = 0 the integral is simplified into 

1 a 

= - e-^x^ + ^+^Ll\{x)dx. (116) 

2 Jo 


Substituting these results in Eq. (100), for the case 
with Ui > 0 and rif = 0 we find 


Using Eq. (104) with c=l, A = ai, a — 1 = ^-|-1^-|-1, 


i.e. a = ^ ^ + 2, and applying the definition (105) 


j{E2) ^ 


2 ' 2 
the result is 


r (f + ¥ + 1) r f “ 5) /.Qo'i n) ii'(^»+ 2 2 i)p7g» I I 


^yriiir^rii + ai + l)r(ay -|- l)r ^ ~ 1 ) 

while in the simplest case of = 0 and nf = 0 one has 

j{E2) ^ 


r(^ + ^ + 2). 


(117) 


r(f+ ^ 


p v'r(a, + l)r(a/ + l)' 


(109) 


In the simplest case of = 0 and n/ = 0 one has 

/(^i)(0,0) = h(| + ^+2). (118) 


In the case of rij = 0, n/ > 0, following the same steps 
one finds 


Substituting these results in Eq. (113), for the case 
with rij > 0 and nj = 0 we find 


/(^2)(0,n/) = 

i r(^/ + ¥-f-i) - 
2 n/!r(f^-f-i) ' 


a/ 

2 


( 110 ) 


j(.El) ^ 
d 


r(f+ ¥ + 2)r(n, + f--^-1) 


(119) 


j{E2) ^ 


y/nP.Tin,+a, + l)T{af + l)r 
while in the simplest case of = 0 and nf = 0 one has 


r(f+ ¥ + i)r("/ + ¥-^-5) 

^n/!r(n/ + af + l)r(ai -I- i)r {°y ~ ^ ~ h) 


( 111 ) 


j{El) ^ 1" (^ + ¥ + 2 ) 


p v/r(a, + l)r(a/ + l)' 


( 120 ) 
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In the case of = 0, n/ > 0, following the same steps 
one finds 




1 r (rt/ + ^ - f - 1) ^ 

2 n;!r(^-f-l) V2 2 J’ 

( 121 ) 


r(Bl) ^ 
P 


r (f + ^ + 2) r (n^ + f - f -1) 

v/n/!r(n/ + a/ + l)r(a, + l)r (1^ - ^ - l) ' 


( 122 ) 


C4. Normalization 

The total wave functions are given in Eqs. (§ and 
(20). The integration over the Euler angles 9 and the 


relevant normalization have been studied in Appendix 
B, while the rest of the integrations are performed over 
/ / /32'^/32/33d/33, where the Z?!, Z?! factors come from the 
volume element and cancel with the first factor of Eq. 
( [^ . Using Eqs. @ and ([^, as well as the relevant Ja¬ 
cobian, one finds that the rest of the integrations are over 
f f fi/32d/33 = J f Pdpdfj). The integration over (j) 

and the relevant normalization factors have been stud¬ 
ied in Appendix A. We determine here the normalization 
factors related to the (3 integration. We have 


1 _ 2(n)! B 

A^l r(n -|- a -|- 1) B2B3 


I{n,n) (123) 


with 


/(n,n)= / Ll0^)Ll0^)d'^. (124) 

Jo 

Using the substitution 0 = x with dx = 2pd$, the inte¬ 
gral is written as 

I{n,n) = ^J e~^x°-Ll{x)Ll{x)dx. (125) 

Considering the case with n = 0, which is of interest 
here, this integral is of the form of Eq. (1041 with c = 1, 
a — 1 = a, thus leading to 


7(0,0) = 


r(a- 7 l) 


(126) 


Then Eq. (123) for n = 0 leads to 
1 B 




(127) 


Appendix D. 


Dl. Kinetic energy and volume elements 

The expressions for the kinetic energy and the volume 
element depend on the dimensionality of the space con¬ 
sidered. We distinguish three cases, with dimensionality 
five, four, and three respectively. 

1) In the usual Bohr Hamiltonian describing the 
quadrupole degree of freedom in the five-dimensional 
(5D) space of the collective variables /3 and 7 and the 
three Euler angles {9, (f), ip), the kinetic energy term reads 

m 


,^)vih 

9? \1 d ^ d 1 d . 80 

2B dp dp /32 sin 37 97 '^d^^ 


(128) 


resulting from the Pauli-Podolsky quantization proce¬ 
dure m in the full 5D space. The volume element reads 

m 

dr = P'^l sindql sin0 dP dj d9 dp dip. (129) 


If the P variable is separated from the rest, either exactly, 
as in the E(5) critical point symmetry |I8) . or through 
an adiabatic approximation, as in the X(5) approach [I6j . 
the volume element in the P part of the problem becomes 

[iniiiH] 

dT0P)=P^dp. (130) 


2) In the Davydov-Chaban approach [7T] the 7 vari¬ 
able is removed from the Hamiltonian from the very be¬ 
ginning of the problem, since 7 is treated as an effective 
deformation parameter. Then the quantization proce¬ 
dure is applied in the 4D curvilinear space of P and the 
three Euler angles. As a result the kinetic energy term 
of the Hamiltonian is obtained in the form 


- {^)vih 


fp_ 

2B 


J_A/33A 

P^dp''^ dp 


(131) 


Note that now the power of /3 in (131) is 3 and not 4 as in 


the /3-part of (128), while the respective volume element 


IS 


dT(^p)vib{P) = P^dp. (132) 

If the wave function is sought in the form m 

iPiP) = P-^I\{P), (133) 

the kinetic energy term in the Schrodinger equation for 
the wave function 0P) appears in the form |7T] 


This result indicates that when calculating /3-integrals in 
B{EL)s, the normalization factors cancel out with 
the g factor appearing in the volume element and 
therefore do not affect the final results. 


~ _ 92 ^ -itP 

(/ 3 )mb - + 8H2/3I ’ 


(134) 


where the second term is further moved into the effective 
potential part. 
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3) In the limit of strong 7 instability of the Wilets-Jean 
approach [72] the nucleus is considered as a droplet which 
can only execute axially symmetric vibrations. This sys¬ 
tem has only three degrees of freedom: /3, 0 and (j). Then 
the kinetic energy term in the Hamiltonian becomes 


D2. Moments of inertia 

Using the standard Bohr expression for the nuclear sur¬ 
face in the body-fixed frame, given by |69| 




(135) 


R{9, = i?o 


Im 


(141) 


where wave functions of the form ipiP) = P 
considered and the volume element reads 


cLt = /3^sin0 dp dO dp. (136) 


This approach has been recently used in Ref. m- 
The kinetic energy term of the Davydov-Chaban ap¬ 
proach has been generalized from quadrupole to any mul¬ 
tipolarity A by Williams and Davidson m, the final re¬ 
sult being 


Tx 


11 5 3 a ■ 

^ [MWxWx. ■ 


(137) 


The basic assumption behind this derivation is the re¬ 
quirement of no vibration-rotation cross terms OESI, 
which diagonalizes the inertial tensor and hence the ro¬ 
tated coordinate system is the principal inertial (body) 
system. 

In the case of simultaneous presence of quadrupole and 
octupole deformation, the kinetic energy within this gen¬ 
eralized Davydov-Chaban approach reads 




E 


A= 2,3 


tP 1 9 3 9 

^xWxWx^^Wx' 


(138) 


Using wave functions of the form 

P^{P2.P^) = {P2P^)-^'\^{P2.PZ), 

which is a straightforward generalization of Eq. 
the kinetic energy takes the form 


(139) 

ra. 


f - - _i_ 

^^^^\2BxdPl^ SBxPl)' 


(140) 


where again the second term is pushed into the effective 
potential, as in Eq. (134). 

Erom the considerations given above, it becomes clear 
that the kinetic energy term used in Refs. [HUSH ED, 
as well as in the present work, is based on the following 
assumptions: 

1) The 7 degree of freedom is frozen from the very be¬ 
ginning, thus reducing the degrees of freedom to four (P, 
three Euler angles) in the case of pure quadrupole de¬ 
formation, and to five [P 2 , Ps, three Euler angles) in the 
case of simultaneous presence of quadrupole and octupole 
deformations. 

2) Vibration-rotation cross terms are ignored, making 
the inertial tensor diagonal and allowing the rotated coor¬ 
dinate system to be the principal inertial (body) system. 


where Yim{ 0 ,p) stands for the spherical harmonics, ig¬ 
noring vibration-rotation cross terms as above, and as¬ 
suming that only the even components ( 030 , a 3 ± 2 ) of the 
octupole parameters are non-vanishing, we obtain for the 
moments of inertia in the octupole degree of freedom the 
expressions [73 [HI 

= 53(6030-1-27^^030032 + 8032) ( 142 ) 

^2^^ = 53(6030 — 2-^/^030032 + 8032) (143) 

= 8530^2, (144) 

which in the axial case (030 = P 3 , 032 = 0 ) give 

af) = = GB3PI = o. (i45) 


Different expressions for the moments of inertia are ob¬ 
tained if one considers the odd components (o 3 ±i, 03 ^ 3 ) 
as the non-vanishing ones m- Here we make the as¬ 
sumption, as in Ref. [73176], that for low-lying collective 
negative parity states the even components play the main 
role, since their contributions to the shape are more sym¬ 
metric, a property usually associated with lower energy 
configurations. 

For the moments of inertia in the quadrupole degree 
of freedom we use the standard expression [53 


df ^ = 452/32 sin^ ( 7 - -nk 


(146) 


which in the axial case (7 = 0) gives 


■vli —32 “ 3i)2P2, 


+ 2 ) 


= 0 . 


(147) 


Collecting (145) and (147) into the axial quadrupole- 
octupole moment of inertia one gets 


( 2 -H 3 ) ^ ^ -j( 2 + 3 ) ^ 




= 3B2Pi + GB3Pi, (148) 


which gives exactly the denominator in the angular mo¬ 
mentum part of Hamiltonian 0 

2dW ^ 6(52/31 + 253/3|) ' 

From the considerations given above, it becomes clear 
that the moment of inertia term used in Refs. [HEa, 
as well as in the present work, is based on the following 
assumptions: 

1) Vibration-rotation cross terms are ignored, as in the 
case of the kinetic energy. 
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2) Only the axial components of deformation are taken 
into account, both in the quadrupole and in the octupole 
degree of freedom, based on the qualitative expectation 
that more symmetric configurations would lie lower in 
energy. 

It should be noted that in Ref. m an expression 
3(732/31 + 733/3|) has been used for the moment of inertia. 


D3. Separation of variables 


Exact separation of the /3 and 7 variables in the frame¬ 
work of the Bohr Hamiltonian can be achieved by con¬ 
sidering potentials of the form u(/3,7) = u(/3) -I- w(7)//3^ 
[711177]. In contrast, when the potential is of the form 
ui/3, 7) = v{j3) + ui(7), only approximate adiabatic sepa¬ 
ration of variables can be tried, as in the case of the X(5) 
critical point symmetry [161 [78]. In the case of X(5), a /3^ 
term survives in the differential equation involving the 7 
variable, replaced in the adiabatic approximation by the 
average value (/3^). The accuracy of this approximation 
has been tested in Ref. m and the limits of its validity 
have been pointed out. The recently developed Algebraic 
Collective Model [HDH55] offers a path for avoiding this 
approximation by performing rapidly converging exact 
numerical calculations instead of pursuing approximate 
analytical solutions. 

In the present case, the 7 variable has been “frozen” 
from the very beginning, following the Davydov-Chaban 
approach m, as explained in Appendix Dl. Therefore, 
no question of separating the /3 and 7 variables appears. 
However, separation of the /3 and cj) variables is desirable, 
in order to achieve analytical solutions in closed form. By 
analogy to the X(5) situation described above, a poten¬ 
tial of the form v{j3,(j)) = u{(}) -I- w{ip^) has been cho¬ 
sen and adiabatic separation of variables has been tried, 
taking advantage of the fact that the w{(l)^) potential is 
supposed to be of the form of two very steep harmonic 
oscillators centered at the values ±0o- Because of the 
steepness of the oscillators it is plausible to use the adia¬ 
batic approximation in the differential equation involving 
the /3 variable [Eq. ^], by replacing the variable (j) by 
±00- Again in analogy to the X(5) case mentioned above, 
a (/3^) term remains in the differential equation involving 
the (j> degree of freedom [Eq. ([^]. There is no need to 
explicitly determine (/3^), since it enters the parameter b 
[Eq. (15l], determined from E2 transitions as described 
in subsection HI.B.I. 

In other words, we exploit for the separation of vari¬ 
ables the fact that the potential is supposed to 

be of the form of two very steep harmonic oscillators 
centered at the values ±0o- This makes the adiabatic 
approximation of (j) by zL^iq plausible, isolating the two 
very steep harmonic oscillators in the (j) equation and 


leaving the rest of the terms in the P equation. An al¬ 
ternative possibility is to consider a potential of the form 
v{P,(j)) = u{p) -\-'w{p^)/P'^. Then the separation of vari¬ 
ables will become exact, but the distribution of terms in 
the two equations will be different. 

The adiabatic approximation used here, based on two 
very steep harmonic oscillators, does have a cost. It is 
well known that the correct description of the parity split¬ 
ting, usually depicted as the odd-even staggering of the 
energy levels of the ground state band and the negative 
parity band, requires a finite barrier between the two 
wells, which is angular momentum dependent |55H57) . 
In the present approach a practically infinite barrier be¬ 
tween the two wells is used for all angular momenta. This 
has as a consequence that the theoretical predictions for 
the low lying negative parity states (especially for 1“ and 
3“) are poor, as pointed out in subsection HI.A. 

It should be mentioned that the Bohr Hamiltonian has 
been solved for the potential 1/ sin^(37) [resembling the 
last fractional term in Eq. 0]. possessing a minimum at 
7 = 7r/6, first by replacing the 7 variable in the moments 
of inertia by its expectation value, 70, and subsequently 
avoiding this approximation |84j , the results revealing the 
approximation to be a good one. Future tests of similar 
nature in the present framework are desirable. 


D4. Comparison to other approaches 


As it has already been mentioned in the Introduction, a 
more general approach has been developed by Bizzeti and 
Bizzeti-Sona m HO], in which nonaxial contributions, 
small but not frozen to zero, are taken into account. It 
is worth commenting briefly on the relation between the 
two approaches. 

1) In the AQOA approach, no nonaxial contributions 
are taken into account. As a result, all variables related 
to nonaxiality in Ref. |19j are vanishing, the matrix of 
inertia (Table I and H in Ref. [19]) becoming diagonal. 

2) Because of the same reason, in the invariants up to 
fourth order reported in Table VHI of Ref. m, only the 
first term in each invariant, containing only P 2 and/or P 3 , 
is surviving. In the present approach only the invariants 
up to second order, being equal to /3| and /3|, are used. 

3) In Ref. [19], in addition to the infinite square well 

potential, a harmonic oscillator potential proportional to 
the square of a; = P 2 has been used. In the present 

approach, P^lP2 = ta^4>y^B2/B3. Therefore the two po¬ 
tentials coincide, up to constant factors, for small angles, 
for which tan p « sin p k, p. 

4) The total moment of inertia appearing in Ref. [TS] 
[Eq. (32)], coincides with the total moment of inertia 
used here [Eq. (149)], if the nonaxial variables vanish, as 
seen from Eq. (6a) of Ref. [T^] . 
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